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Abstract 

This paper covers the development of an integrated nonlinear dynamic simulation for 
a variable cycle turbofan engine and nozzle that can be integrated with an overall vehicle 
Aero-Propulso-Servo-Elastic (APSE) model. A previously developed variable cycle turbo- 
fan engine model is used for this study and is enhanced here to include variable guide vanes 
allowing for operation across the supersonic flight regime. The primary focus of this study 
is to improve the fldelity of the model’s thrust response by replacing the simple choked flow 
equation convergent-divergent nozzle model with a MacCormack method based quasi-lD 
model. The dynamic response of the nozzle model using the MacCormack method is veri- 
fled by comparing it against a model of the nozzle using the conservation element/solution 
element method. A methodology is also presented for the integration of the MacCormack 
nozzle model with the variable cycle engine. 
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Nomenclature 


Cross-sectional area 
Aero-propulso-servo-elasticity 
Aero-servo-elasticity 
Conservation element /solution element 
Convergent divergent nozzle 
Computational fluid dynamics 
Courant-Eriedrichs-Lewy stability condition 
Computational fluid dynamics flux vector 
Eully unstructured Navier- Stokes in 
three dimensions 
Source term vector 
Length 

Numeric propulsion system simulation 

Mach number 

Mass flow parameter 

Pressure 
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Temperature 

Solution vector 

Volume 

Variable cycle engine 


VGV Variable guide vane 

W Mass flow 

a Speed of sound 

Cy Specific heat at constant volume 

e Internal energy 

t Time 

u Velocity 

X Length 

7 Ratio of specific heat 

p Density 

Subscripts 

N Total number of grid points 

be Boundary condition 

i Station location 

s Static flow condition 

t Total flow condition 

Superscripts 

' Non-dimensional variable 

* Throat location 
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I. Introduction 


T he NASA Fundamental Aeronautics Program, High Speed Project, aims to overcome the obstacles 
associated with supersonic commercial flight. This project aims to develop the technologies to allow for 
practical supersonic overland commercial transport at a Mach number slightly less than two. The proposed 
vehicles are long, slim body aircraft with the potential for pronounced structural vibrations known as aero- 
servo-elastic (ASE) modes. When coupled with propulsion system dynamics, the modes are known as 
aero-propulso-servo-elasticity (APSE). The APSE considerations can lead to design challenges pertaining to 
aircraft performance such as aircraft ride quality and stability. Eurthermore, other disturbances upstream 
of the inlet generated by atmospheric wind gusts may also affect ride quality and stability. To study the 
relevant dynamics of the propulsion component of APSE, a nonlinear variable cycle engine (VCE) model 
was developed^ and is improved upon by incorporating a computational fluid dynamics (CED) model of the 
convergent- divergent (CD) nozzle. 

CED models of a supersonic CD nozzle are being developed and integrated into a VCE model to capture 
relevant thrust dynamics. Accurate representation of the thrust dynamics is critical for this research effort, 
as these models will ultimately be integrated into an aero-elastic vehicle dynamic model to investigate the 
physical coupling of the integrated vehicle and closed loop control performance. 

This paper investigates the CED modeling of a CD nozzle for a supersonic commercial transport using 
the MacCormack method. Comparisons of the steady state and dynamic result of the MacCormack method 
are verified against a Conservation Element /Solution Element (CESE) CED approach for capturing the key 
dynamics of the CD nozzle. Integration of the MacCormack nozzle model into the previously developed 
VCE model is done to improve upon the simple choked flow equation that was previously used to model 
the nozzle. This higher fidelity modeling of the exit flow of the engine will allow for accurate capturing 
of the thrust dynamics required for the APSE task. The modeling approach used for the VCE model has 
previously been developed^’ ^ and verified.^ The work presented here for the VCE, however; expands upon 
the previous single design point to other key portions within the supersonic operating envelope through the 
implementation of variable stator vanes and exit nozzle areas. The work presented in this paper is only a 
portion of the broader dynamic modeling of the propulsion system for the needs of APSE.^ 

The overall paper objectives are as follows: 

• Provide an overview of the supersonic CD nozzle model, with a focus on integration of the nozzle with 
the VCE propulsion system. 

• Use of the CESE method for comparing dynamic modeling accuracy of the CD nozzle model. 

• Illustrate the key features of the propulsion system model required to accurately capture the thrust 
dynamics. 

The paper is organized as follows. Eirst, a modeling overview that includes descriptions of the overall 
APSE goals, the CD nozzle, and the VCE model with a focus on variable guide vanes and nozzle integration. 
Second, a results section that will focus on the dynamic accuracy of the nozzle model and its integration 
with the VCE model for accurate thrust estimation. Third, the paper will outline the next steps required 
to complete the overall dynamic propulsion system, primarily the modeling of the external compression 
inlet, and its integration into the overall elastic vehicle model being developed in NASA Langley’s fully 
unstructured Navier-Stokes in three dimensions (EUN3D) code. Einally, some concluding remarks will be 
provided on the performance of the integrated VCE- nozzle model. 

II. Modeling Overview 

A. Supersonic APSE Overview 

The overall APSE simulation block diagram can be seen in Eig. 1. The propulsion system of the overall 
model, highlighted in Eig. 1, is comprised of the inlet, VCE, CD nozzle, and propulsion controllers. ASE 
vehicle impacts are neglected in this study, but are illustrated by dashed lines in Eig. 1 to provide an overview 
of how the propulsion system fits into the larger APSE simulation architecture. A mixed compression inlet 
was used in previous studies.^ However, the current project design is for an external compression inlet and 
will be modeled in this study using simple steady state performance calculations. All of the other propulsion 
elements are modeled using unsteady conservation equations, described in the following sections. This allows 
for the investigation of VCE control performance and thrust oscillations. To run the simulation, free stream 
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Flow Perturbations 
from Vehicle 
Vibrations 


static conditions of pressure, temperature, and Mach number are applied to obtain an inlet exit total pressure 
and temperature, which are then used as the engine face input conditions. The engine then uses the input 
conditions, its component performance maps, and a modified version of the Euler conservation equations 
to calculate the engine output provided to the CD nozzle interfaces. The version of the VCE used for the 
current study contains two CD nozzles, one for the main flow through the core and primary bypass called the 
primary duct nozzle, and a second for the variable cycle duct exit called the VCE nozzle. Both the primary 
duct and VCE duct nozzles are modeled using the MacCormack method. 

“V As discussed, an overall goal 

of the ASE project is to integrate 
a high fidelity vehicle and propul- 
sion system model to investigate the 
dynamic performance coupled with 
the elastic vehicle. The work pre- 
sented here is a step along that 
path, providing improved dynamic 
thrust estimation. Previous work in 
the area of APSE has focused on ei- 
ther the vehicle or the propulsion 
system, but rarely are both systems 
modeled with the same level of fi- 
delity. In addition, most of the pre- 
vious research has only considered a hypersonic vehicle where the vehicle itself provides compression of 
the flow feeding into the propulsion system. This leaves the speed regime of supersonics and the use of 
turbo-machinery as the propulsion system as an open field of research. 

B. Nozzle Model 
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Figure 1. Block diagram of overall APSE system model with propulsion 
system highlighted. 


A CD nozzle was modeled to investigate the dynamic thrust behavior caused by upstream flow perturbations. 
As previously mentioned, the VCE propulsion system has two nozzles, one for the exit of the core and bypass 
duct called the primary duct nozzle and one for the variable cycle diet exit called the VCE duct nozzle. Both 
nozzles are modeled exactly the same with only geometric differences. The primary nozzle chosen for the 
supersonic transport includes a variable geometry CD internal duct and a conical plug or center body. The 
CED approach for this study is geared towards having as simple of an approach as possible while still 
capturing the relevant thrust dynamics, thus a quasi- ID method is used. Thus, the approach described here 
is only suitable to model the internal duct portion of the two nozzles and neglects any effects of the external 
plug. This approach is justified by preliminary studies using a 2D CESE code that showed the external 
portion of the nozzle only contributed a few percent to the overall thrust. In addition, the ultimate goal of 
this work is to have a propulsion system model that will be integrated into NASA Langley’s EUN3D^^ code. 
This is a widely used 3D CED code, where the propulsion system presented here will provide the dynamics 
of the internal duct and EUN3D will handle the more complex external flows. 

Both a CESE and MacCormack quasi-lD method are used for modeling the internal duct of each nozzle, 
given the following Euler equations for mass, momentum, and energy conservation respectively: 


d(pA) du dA dp 


du du ^ 


^dx dx 


dT dT r du dln(A) 


( 1 ) 

( 2 ) 

( 3 ) 


The above equations Eq. 1-3 are perfectly fine to implement into a numerical scheme. However, for this 
study a non-dimensional form of the equations will be used. This is more style than substance, but can help 
to visualize changes in parameters without having large scaling issues for visualization. The dimensionless 
terms chosen are described below and denoted with a prime, 
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CH 
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A* 
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The conservative form of the governing equations was chosen for implementation. There is no mathemati- 
cal difference from the non-conservative form, however the use of the conservative form can provide improved 
numerical results when the flow has discontinuities, such as the case of internal shocks in the nozzle. The 
conservative form can be written in a compact vector form as: 


dlJ _ dF 
dt' dx' 

where the components of the solution vector (U), fluxes (F), and source (J) are: 

Ui = p'A', U2 = p'A'u', Us = p'A 


^ 7 /2 


Fi = p'A'u', F2 = p'A'u'‘^ F -p'A', Fs = p'u'A' 


7 


7-12 


7-1 2 

p' A'u‘ 




1 BA' 

Ji=0, = Js = 0 

7 Bx' 


(10) 


( 11 ) 

(12) 

(13) 


The two modeling approaches allow for comparison and simple verification purposes. For the purposes 
of this paper the MacCormack method will be used as the primary modeling approach, meaning it will be 
the method that is integrated with the VCE model. The CESE method is an independent code used in this 
study for comparisons with the MacCormack method to improve the level of confidence of the nozzle model. 
The MacCormack method was chosen for integration with the VCE model for its simple finite difference 
approach. This approach more closely resembles the modeling method utilized for the engine. The following 
subsections will provide more details of the two approaches. 


MacCormack Method 

The MacCormack method is an explicit finite difference scheme that is very popular in CED. This scheme 
allows for the accuracy of a Lax-Wendroff scheme, while reducing the overall complexity. The MacCormack 
method is multi-level, meaning that there is a multistep process to reach the solution. These steps are 
traditionally called predictor and corrector, where an averaging can then be performed. The averaging at 
each time step is what allows for the higher accuracy while using a relatively simple finite difference scheme. 
A step-by-step procedure of the MacCormack method can be found in Anderson. 

In the subsonic portion of the nozzle, one of the boundary condition variables is required to float, since 
information under this condition traverse both upstream and downstream. The choice was to allow the 
velocity to float to be consistent with the method outlined in Anderson. The boundary conditions are then 
defined in terms of the solution vector, where the density and temperature at the nozzle entrance are fixed: 


^1,2=1 = PbcAi 

(14) 

U2,i=l = 21/2,2=2 — 1/2, 2=3 

(15) 

TT ( 1 7 2^ 

"=^-^77-1 2^ j 

(16) 
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The boundary conditions at the nozzle exit are extrapolated using a simple linear relation, since the flow 
is assumed to have no internal normal shocks, and thus be supersonic: 


U\^N = — U\^n-2 (17) 

^2, AT = 2[/2,AT-1 — U2,N-2 (18) 

Us,N = 2t/3,AT-l — Us^N-2 (19) 

Computational Element / Solution Element Method 

The CESE method is a second-order accurate method that has been applied to a wide variety of CED 
problems and has been proven to be a simple, robust, and efficient tool for dynamic simulations and aero- 
acoustic computations. A guiding principal of the CESE formulation is that space and time are unified 
and time is simply treated as an additional dimension to the problem. In addition, the enforcement of flux 
conservation in space and time at an interface is an integral part of the solution procedure. The method also 
uses a staggered mesh. The CESE method applied here for the modeling of a quasi-lD CD nozzle, expands 
upon the work by Chang for the modeling of a shock tube.^^ The approach here follows Changs developing 
of the method applied to the Euler equations with the primary difference being that the equations above 

are used for a CD nozzle instead of a constant area shock tube as presented by Chang. This leads to the 

governing equations in the previous work needing to be modified to account for the source term that exists 
due to the area variation. Otherwise the exact same procedure is followed. 

The boundary conditions for the CESE method are applied as outlined above in the MacCormack method, 
by simulating a subsonic entrance to the nozzle and supersonic exit. A slight variation to the application of 
the boundary conditions is that the interpolation and extrapolation need to account for the staggered grid. 

C. Engine Model 

The current engine de- 
sign concept for the super- 
sonic propulsion system is a 
variable cycle engine (VCE), 
which is an extension of a tur- 
bofan engine. The concept 
behind a VCE is to optimize 
the cycle throughout the flight 
regime by trimming the flow 
through the VCE fan stream 
so that the overall engine flow 
demand matches that of the 
inlet capture air. This is done 
in order to maximize pressure 
recovery and minimize flow 
spillage, which adds to drag. 

The VCE fan stream is also utilized for noise abatement for supersonic overland flight. 

The flow path of the VCE is split into three gas paths as illustrated in Eig. 2. The primary gas path is 
nearly identical to a turbojet and goes through the core of the engine. The secondary gas path is similar to 
a typical turbofan engine in that a large amount of the flow is bypassed around the core of the engine. In 
the VCE version covered in this paper, the main bypass flow is mixed with the core flow downstream of the 
turbine and exhausted out the common primary duct nozzle. The more advanced concept part of this engine 
design is a third gas path that exhausts through the VCE duct nozzle. The fan in Eig. 2 is split, such that the 
flow going to the bypass and core is one component model called fan, and the fan tip region feeding the third 
flow path is modeled by a component call VCE fan. The third flow path provides a lower exit velocity exiting 
the VCE duct nozzle that could be used as a noise shield for the flow exiting the primary duct nozzle. While 
Eig. 2 shows all of the major engine components that have their own conservation equations, in between 
each component is also an associated connecting duct that is not shown in the figure for simplicity. Each 
of these ducts are also modeled using the governing conservation equations.^ In addition, a preliminary fan 
speed controller is included that was previously developed. 
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Figure 2. Engine design concept for a variable cycle engine. 
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The initial operating conditions and component geometry information, as well as the component per- 
formance maps required to implement the dynamic VCE model is obtained from a design concept of a 
supersonic VCE implemented in the Numerical Propulsion System Simulation (NPSS).^^ A cruise operating 
point was chosen at 15,240 m (50,000 ft) and a Mach number of 1.6. A choked nozzle flow equation has 
previously been used for the exit boundary condition. The VCE model uses a single lumped volume for each 
of the major components such as compressor, combustor, turbine, afterburner, and nozzle. The modeling 
approach is outlined in previous works. ^ A, is Each of the fluid flow components is modeled using a set of 
derived conservation equations modified from the standard Euler form and written for continuity, momen- 
tum, and energy. These equations are integrated numerically using a time marching scheme and the Seldner 
differencing technique of the spacial terms. ^ 

In typical subsonic commercial aircraft, CED models of turbofan engines have not been necessary for 
control design due to the operating conditions and expected lower frequency perturbations. Eor such models 
in this operating regime, only the turbofan shaft inertias have been required for closed-loop controls devel- 
opment. However, for new commercial supersonic transports the assumptions of just being able to develop 
propulsion system controls around shaft dynamics may no longer be valid. 

Variable Cycle Engine Variable Guide Vane Modeling for Expanded Elight Envelope 

The primary update of the dynamic VCE model for the work presented here from previous studies is the 
inclusion of variable guide vanes (VGV) that allow for the engine to operate across its full supersonic operating 
envelope. Previous models where designed about a single operating point. However, numerical instabilities 
could be encountered when attempting to operate these models at other points throughout the supersonic 
operating envelope. The VGV model is developed and integrated with the compressor components, which 
include the high pressure compressor, fan, and VCE fan. This allows for variable geometries necessary to 
operate the VCE within its operating range of altitude, Mach, and power, while forcing the compressors (via 
control schedules) to operate on a desired operating line of the compressor performance maps. 

The VGV model consists of performance characteristics with sets of maps at different VGV angles, 
which are interpolated. An example of the VCE fan performance map and its variation with VGV angle 
is illustrated in Eig. 3. The solid black lines represent corrected speed and the color map represents the 
operating efficiency. Eor this dynamic model, the effective cross sectional area is calculated as a function 
VGV angle. The maps supplied by the NPSS model include the overall incidence angle. The incidence 
angle, however, is taken out in this model and the reference angle with no VGV turning is referred to as 
the nominal angle. This nominal angle corresponds to the initialized compressive component cross sectional 
area. Eor the VGV model, the mass flow parameter {Mpp) is utilized and the compressor effective flow area 
is calculated as a function of the Mpp. Additionally, by taking into account the cosine of the VGV angle as 


MFpPtCos(a) 

Mfp = 

\ R 2 

where A, Mpp^ kE, Pt^ Tt, cos((a), and M are the component cross section area, the mass flow parameter, the 
component total pressure and temperature, the cosine of the VGV angle, and the component Mach number. 
The angle, a, is the VGV turning angle. 


( 20 ) 

( 21 ) 


NAS A/TM— 20 15-21 8479 


6 



Figure 3. The VCE fan compressor performance maps at a low VGV angle (left) and a high VGV angle (right), 
where the color map corresponds to efficiency. 


D. Engine and Nozzle Integration 

The integrated nonlinear propulsion system used in this study includes the VCE and CD nozzle models. 
This allows for the development of the modeling approach, control concepts, and interfaces for the overall 
propulsion system integration and APSE interfaces. Ultimately, an external compression inlet needs to be 
developed and integrated with the VCE to have a more accurate picture of how atmospheric disturbances 
will impact thrust for the current supersonic design concepts. This effort is currently ongoing.^ 

The integration approach is to replace the previous choked flow equation of the engine model exit bound- 
ary condition with the MacCormack method CED CD nozzle described in the preceding sections. The choked 
flow boundary condition is replaced with a very small duct volume that can be added to the engine simula- 
tion. This duct volume is modeled utilizing volume dynamics, similar to the approach applied in modeling 
the ducts in the VCE, where the conservation equations become Eqs. (22) to (24). 

_ ^nozzle^^ ( 22 ) 

-^(to) = _ pvccM^ (23) 

j^iPsTt) = ( 24 ) 

This duct volume serves as the VCE- nozzle interface used to calculate the states and the static density 
and temperature required for the nozzle model. The MacCormack nozzle then provides the mass flow rate 
that was previously supplied by the choked flow equation. The superscripts used in the above equations 
indicate the location of the variable relative to the VCE-nozzle interface. These equations are then tied 
together using the state equation. The nozzle model then provides the dynamic estimation of gross thrust. 

E. Engine Fan Speed Controller 

The engine control design used in this study has aspects of classical loop shaping^^ design and quantitative 
feedback theory. The engine controller for the fan speed was first designed based on linear models, and 
later implemented in the nonlinear VCE model. The fan speed is controlled using the engine fuel injector 
that has a bandwidth of 6 Hz. In addition, the VCE engine design has variable nozzle and VGV that are 
set on schedules and not actively controlled. The nozzle geometry and VGV is currently scheduled based on 
the corrected core speed. 
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III. Results 


The purpose of these simulation results is to demonstrate the steady-state accuracy and dynamic perfor- 
mance of the CD nozzle model using the MacCormack method and the integrated propulsion system model 
with a fan speed controller. The results will be broken out into two sub-sections; one for the stand-alone 
nozzle component verification and the other for the nozzle integrated with the VCE. Each of the sub-sections 
will first discuss the steady state accuracy and then focus on the dynamic performance. The VCE-nozzle 
integration sub- section will focus on the thrust dynamics. 

A. Nozzle 

This sub-section will investigate the steady state and dynamic results for the primary duct nozzle. The 
VCE duct nozzle model is also included for the final integration with the VCE, however, for brevity only 
the primary duct results are presented here. Comparison with the CESE method are provided and further 
studies are done with the MacCormack method, since this is the method utilized for the final nozzle model 
that is integrated with the VCE model. The results presented here are for a cruise condition of Mach 1.6 
and 50,000ft. 

Steady State 

The first comparison of the Mac- 
Cormack and CESE nozzle are sim- 
ply the steady state results com- 
pared to the known analytical solu- 
tion based on the area-Mach num- 
ber relationship.^^ Steady state re- 
sults from these simulations for the 
primary duct are shown in Eig. 4. 

The comparison is done across the 
length of the nozzle for the non- 
dimensional primitive variables of 
density, velocity, and temperature. 

The MacCormack method is shown 
as solid lines, the CESE method 
as dashed lines, and the analyti- 
cal solution as “x’s”. It can be 
seen that the CESE solution is in 
excellent agreement with the Mac- 
Cormack method and analytical so- 
lution. The only slight deviation 
shown is near the throat of the noz- 
zle, located at an axial location of 
0.55. 

A r , 1 ,1 . j-v, i\/r Figure 4 . Nozzle steady state comparison for MacCormack, CESE, and 

IU.I*ljR.GI* S uU-CX\^ U-Sm£f uR-G l.vJ.R)C“ ij_* 1 ij_* I? 1* 1 • I • 1 . . | 

^ analytical solution tor normalized static density, temperature, and ve- 

Cormack method was conducted to locity. 
determine the steady state accuracy 
with various grid sizes compared to 

the analytical solution. The results are tabulated in Table. 1. The internal grid points, neglecting the 
boundary conditions, are used to calculate the mean percent error across the entire computational domain 
of the nozzle. The steady state results show that even with a very sparse grid of 21 points, the mean error is 
less than 1%, and substantially decreases with larger grid sizes. Note, while the grid sizes where increased, 
the Courant-Eriedrick-Lewy (CEL) stability condition was held constant at 0.95, which essentially changes 
the time step of the simulation as well. 
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Table 1. Nozzle steady state mean absolute percent error with various grid size using the MacCormack method. 


Primitive Variable 

21 pts 

101 pts 

201 pts 

Velocity 

0.10 

0.0088 

0.001 

Temperature 

0.21 

0.0026 

0.00014 

Density 

0.67 

0.016 

0.00289 


Dynamic 

For the dynamic results, a grid study of the MacCormack method was conducted to determine the relative 
number of nodes required to get minimal changes in the magnitude and phase of a bode plot spanning the 
frequencies of interest. The highest expected frequency due to atmospheric turbulence at the desired cruise 
condition of less than Mach 2 is about 20Hz. However, the resulting aero-elastic disturbances due to the 
vehicle vibrations could get as high as 60Hz. This dictates that the frequency range of interest is extended 
to 600 Hz, since dynamics up to this frequency can still effect the phase at 60 Hz. Dynamic results for the 
primary duct using the MacCormack method were generated where various grid sizes are used to investigate 
impact on dynamic accuracy. 

The grid study is done using a sinusoidal logarithmic sweep of 1% amplitude on the nozzle density 
entrance boundary condition, and then compared to the density at the nozzle exit. The resulting bode plot 
is shown in Fig. 5. The left of Fig. 5 shows the full frequency range from 6 to 600Hz, and the right side is a 
close up at the higher frequency range. It is illustrated that even a very small number of nodes can provide 
reasonable accuracy at frequencies of less than lOOHz, but the sparse grids begin to show deviation beyond 
that frequency. From these results it was determined that 101 nodes or grid points would be sufficient for 
the frequency range of interest. This was determined since the variation from the case with 201 nodes at the 
600Hz range was less than IdB in magnitude and only a couple degrees in phase. 
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Figure 5. Nozzle exit density bode plot obtained from a logarithmic sinusoidal density wave applied at the 
nozzle entrance boundary condition for investigating the dynamic accuracy due to grid sensitivity. 


Upon completion of the grid study, the MacCormack and CESE methods were compared at various 
constant frequencies. Similar to the bode plots in the grid study a 1% amplitude disturbance on the nozzle 
density entrance boundary condition is compared to the density at the nozzle exit, but now using a constant 
sine wave. The results of such a comparison at lOOOHz are shown in Eig. 6. The top plot shows the 
normalized disturbance of density subtracted by the steady state value versus non-dimensional time units. 
The bottom plot shows the results of the density response at the nozzle exit for both the MacCormack 
(black) and CESE (red) methods. The results show that the two methods agree very well dynamically even 
at higher frequencies. 
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The constant sine wave disturbance de- 
scribed above is then repeated at various dis- 
crete frequencies spanning the frequency range 
of interest. The input to output response is 
then determined and used to generate a bode 
plot. This is done in order to compare the fre- 
quency response of the MacCormack method 
against the CESE method for discrete frequen- 
cies, with the results shown in Eig. 7. In the 
figure, the discrete results using the MacCor- 
mack method are plotted as blue crosses and 
the CESE method are shown as red circles. 
It can be seen that two methods agree very 
well, and provide a high level of confidence 
in the nozzle dynamics in the absence of ex- 
perimental data. In addition to the discrete 
sine waves, a logarithmic sinusoidal sweep is 
applied using the MacCormack method to ob- 
tain the frequency response shown as a solid 
blue line. The sinusoidal logarithmic sweep is 
a more desirable method to evaluate the dy- 
namic response as it takes much less computa- 
tion time. It is assumed here that the dynamic 
response obtained by the discrete frequencies is more accurate than the sweep. The results shown in Eig. 7 
illustrate that the sweep is able to accurately obtain the magnitude response, with only an error of about 
0.5 dB. The accuracy of the phase is not as good, particularly at the higher frequencies. However, even at 
the higher frequencies the difference between the sweep and constant sine waves is only about ten degrees. 
The most critical frequencies to accurately model are 60Hz and below, so with an error of only ten degrees a 
decade above the frequency of interest this result shows that if a quick evaluation of the dynamics is required 
a sweep can be used with relatively good accuracy. 



-0.5 


Figure 6. Nozzle entrance boundary condition disturbance 
of density using a constant sine wave at lOOOHz normalized 
with the mean removed versus non-dimensional time units. 



Frequency, Hz 


Figure 7. Bode plot of the Nozzle density input to output comparing a 1% amplitude constant sine wave using 
the MacCormack method (blue cross) to the CESE method (red circle). In addition, a logarithmic sinusoidal 
sweep is shown using the MacCormack method as the solid blue line. 
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B. Engine-Nozzle Integration 

The following results are for the VCE model integrated with the MacCormack method model of both the 
primary duct and variable cycle duct nozzles. The following section will outline the steady state comparison 
of the model compared to the initialization data from the NPSS model, the fan speed control response, and 
the gross thrust response. The results presented here are for a cruise condition of Mach 1.6, 50,000ft, and 
100% power, unless otherwise stated. 

Steady State 

The VCE dynamic model is initialized by steady state data provided by the NPSS model of the VCE. The 
steady state percent difference between the NPSS model initialization data and the dynamic model of key 
engine components for mass flow, density, and temperature is illustrated in Table 2. The results show that 
the absolute percent error is typically less than 2%, however the compressive components have slightly larger 
errors. The errors here are most likely due to slight geometry differences between the NPSS and dynamic 
model, and differences in the interpolation schemes for the performance maps. In general the steady state 
results are expected to be reasonable from a controls perspective. 


Table 2. Engine-Nozzle steady state absolute percent error. 


Engine Components 

Mass Plow 

Density 

Temperature 

VCE Fan 

2.65 

4.34 

1.43 

Ean 

0.01 

2.37 

0.70 

HPC 

2.05 

2.61 

4.35 

Combustor 

0.49 

0.51 

0.69 

HPT 

0.94 

0.08 

1.66 

LPT 

0.94 

1.09 

0.99 


Dynamie 


The following results will high- 
light the dynamic response of 
the integrated propulsion system, 
which includes the VCE engine and 
CD nozzles using the MacCormack 
method. The results in this sec- 
tion will again use the previously 
discussed cruise conditions at Mach 
1.6. To begin, the VCE controlled 
response is investigated. The VCE 
has a simple fuel actuator with a 
bandwidth of 6Hz, and fuel flow is 
regulated to control fan rotational 
speed. A more detailed discussion 
of the controller can be found in 
previous papers covering the VCE 

engine development.^ ’^3 ^ 

pie series of commanded fan speed 
steps up and down with a percent 
change of 1% from nominal are used 
to illustrate the controlled response 
and are shown in Eig. 8. The dis- 
turbances in this section are applied 

after 2.5 seconds to allow for the startup transient to settle out. The dashed red line here is the commanded 



Figure 8. Engine-nozzle fan speed commanded steps and corresponding 
fan speed response. 
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fan speed and the solid blue line is the response of the fan speed. This modest control demand is just used to 
illustrate that the propulsion system can operate about and maneuver around the cruise operating condition. 

The primary enhancement to the VCE model from previously published work is the inclusion of the 
VGV for the compressive components that allow for operation about the supersonic flight regime. Operation 
across this flight regime is illustrated in Fig. 9. The plots shown here are again for a commanded fan speed 
up and down of 1% from nominal as was discussed above. The plots are of the resulting thrust response 
normalized by the nominal output thrust for the Mach 1.6 operating condition, and the fuel flow is regulated 
to a specified fan speed. The thrust responses are all very similar, with the only significant difference being 
the steady state thrust output of the engine as expected. In the top plot of Fig. 9, the gross thrust output 
is increased by about 44% to operate at a cruise Mach of 1.2 and a decreased altitude of 44,000 ft. The 
bottom plot shows the reduction of the thrust output in relation to the nominal cruise condition to about 
84% when operating at Mach 2.0 and an increased altitude of 52,000ft. 
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Figure 9. The VCE- nozzle simulated gross thrust response due to steps up and down in fuel flow at Mach 1.2 
and altitude of 44,000 ft (top plot), Mach 1.6 and altitude of 50,000 ft (middle plot) and Mach 1.8 and altitude 
of 52,000 ft (bottom plot). 


The thrust response of the VCE integrated with the MacCormack method for both the VCE and primary 
duct nozzles is now compared with the previously used simple choked flow equation. To begin, the same 
commanded fan speed steps are used at the nominal cruise condition of Mach 1.6 and altitude of 50,000 
ft. This is done to illustrate the lower frequency thrust dynamics and is shown in Fig. 10. The overall 
steps up and down are shown in the left most plot where it can be seen that the responses are very similar. 
The propulsion system with the MacCormack method model of the nozzles has a slightly more significant 
response to the steps up and down in fan speed causing a small steady state offset when using the choked 
flow equation. The right plot of Fig. 10 is simply a close up and shows that the transient responses of the 
initial step up are very similar between the two propulsion system models. 

The thrust response shown in Fig. 10 was not entirely unexpected as the flow in the nozzle is supersonic 
and slower dynamic responses should not show significant differences. However, at higher frequencies the 
nozzle model using the MacCormack method will be able to capture the relevant dynamics that the simple 
choked flow equation cannot. For illustrative purposes, a constant density sine wave perturbation is applied 
at the nozzle entrance at a frequency of 600Hz (the highest frequency of interest) and 1% amplitude. This 
is used to compare the simple choked flow equation to the MacCormack method shown in Fig. 11. Here it 


NAS A/TM— 20 15-21 8479 


12 



Figure 10. Propulsion system normalized thrust response to commanded fan speed steps up and down (left) 
and a close up of the first step up (right). 


can be clearly seen that the perturbation using the choked flow equation flows directly through the nozzle, 
dashed red line, without any delay or added dynamic response. The MacCormack method does account for 
the time delay of the perturbation to flow through the nozzle before it impacts the thrust response and it 
captures some added dynamics of the nozzle. This leads to the two responses being out of phase, however 
the time scale is very small. 

Finally, a pressure perturbation 
is applied at the entrance to the 
engine face using a constant sine 
wave of 60Hz and an amplitude 1% 
of the nominal steady state pres- 
sure. The frequency of the pertur- 
bation was chosen as the highest ex- 
pected aero-elastic disturbance fre- 
quency. Ultimately for the APSE 
studies, the perturbations of inter- 
est will be atmospheric wind gusts 
in the free stream. However for 
the present study, a simple pressure 
perturbation is applied at the en- 
gine entrance face due to the ab- 
sence of an inlet dynamic model. 

Normally the external compression 
inlet should provide some dissipa- 
tion of the free stream perturbation, 
so this could be thought of as more 
severe than applying a 1% pressure 
perturbation to the free stream con- 
ditions. This disturbance is illus- 
trated in the top plot of Fig. 12. 

The bottom plot of Fig. 12 shows the resulting thrust responses of the propulsion system with using either 
the simple choked flow equation or MacCormack method model of the nozzles. The two response are different 
with the MacCormack method having a slightly smaller amplitude response to the initial perturbation, but 
the main take away is that the two response are very similar. While the MacCormack method has been 
shown to capture the relevant dynamics of the nozzle as shown in Fig. 11, when a perturbation is applied at 



Figure 11. Nozzle thrust response due to 600Hz constant density per- 
turbation at the nozzle entrance, comparing the simple choked flow and 
MacCormack modeling approaches. 
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the engine face the dynamics of the engine and the fuel flow controller are dominate. The dynamics of the 
nozzle at this given supersonic condition and geometric design causes the nozzle to not play a major role in 
the thrust dynamics of interest for APSE. An operational region where this could change for the proposed 
flexible supersonic vehicle is the case of a takeoff condition, where the vehicle will be at subsonic speeds and 
the dynamics of the flow through the nozzle could change. 

A simple way to increase the fi- 
delity of the choked flow equation 
model of the nozzle would be to in- 
clude a time delay to approximate 
the flow propagation time across the 
nozzle length. Since the ASE cou- 
pling frequency can be as high as 
60Hz, the dynamic accuracy of the 
phase at this frequency between the 
MacCormack nozzle simulation and 
the choked flow nozzle equation is 
critical for control design. To in- 
vestigate the accuracy of the simple 
time delay compared to the Mac- 
Cormack method, data from Eig. 7 
is duplicated and a straight line ap- 
proximation of the time delay phase 
shift in log scale is shown in Eig. 13. 
Eigure 13 shows that this difference 

in phase at 60Hz is only about 6 de- 
Figure 12. Nozzle thrust response to a constant 60Hz sine wave pres- acceptable as 

sure perturbation oi 1% amplitude at the engine entrance race and the r i i i a 

corresponding thrust response for both the VCE with the choked flow error for closed-loop APSE dy- 
equation and MacCormack method nozzles. namics for control stability. Obvi- 

ously, if the nozzle dynamics can be 

simulated by only utilizing the choked flow equation, this will reduce the complexity and execution time of 
the propulsion system simulation. 
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Figure 13. MacCormack method nozzle compared to a choked flow equation nozzle model with a flow propa- 
gation time delay using discrete frequency density perturbation sinusoids over the frequency range of interest 
to investigate phase change. 
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IV. Future Work 


To complete the overall propulsion system model, a nonlinear dynamic model of the external compression 
inlet is required to capture all of the relevant dynamics that could impact the thrust response. While it has 
been shown here that the nozzle may not play a major role in capturing the relevant dynamic response of 
thrust for the purposes of studies in APSE, it is expected that the inlet will play a significant role. This is 
primarily due to the inlet being the interface of the propulsion system to the free stream perturbations of 
interest in APSE. In addition, the inlet is a large volume where the velocity of the fluid is relatively slow 
compared to that of the nozzle. Moving forward, the MacCormack method for the nozzle will continue to be 
used in the final development of the overall propulsion system to get a full assessment of the thrust dynamics. 
The importance of the nozzle dynamic accuracy could play a role once the inlet is connected to the overall 
propulsion simulation or operation of the propulsion system at flight conditions other than those explored 
in this study. Upon completion of the inlet dynamic model and its integration with the VCE model, the 
overall propulsion system needs to be included in an overall simulation of the aero-elastic vehicle. Work is 
currently underway for incorporating the propulsion system models used in this study into NASA Langley’s 
EUN3D code. 


V. Conclusion 

An integrated engine and nozzle propulsion system component model has been developed with a fan speed 
controller, where the engine components have previously been reported upon. The integrated propulsion 
system model presented in this work is suitable for incorporation into an overall supersonic vehicle aero- 
propulso-servo-elastic (APSE) model where accurate thrust dynamics are of primary concern. A convergent- 
divergent nozzle model has been developed for replacing the relatively simple choked flow equation model 
previously used in the engine model to improve the thrust dynamic response. The nozzle model is developed 
using the MacCormack method and its dynamic accuracy is compared against a nozzle model using a 
conservation element /solution element method. The dynamic comparison of the two approaches showed 
very good agreement. In addition to replacing the simple choked flow nozzle of the engine model, variable 
guide vanes for the compressive components were also included. This allowed for operation of the propulsion 
system across its supersonic flight regime, which was previously not possible. The engine-nozzle integrated 
propulsion system showed that at supersonic cruise when including the higher fidelity MacCormack nozzle 
models perturbations at the engine face do not have a significant impact on the overall thrust response 
compared to the simple choked flow model. Eor the desired cruise operating condition and this specific 
nozzle design, a simple choked flow model for the engine nozzle should be sufficient to capture the relevant 
thrust dynamics for APSE needs, as the dynamics of the other engine components and fan speed controller 
dominate the thrust response. 


References 

^Cohen, P., Long-Davis, M., and Povinelli, L., “Fundamental Aeronautics Program - Supersonics Project,” 
http:/ /www. aeronautics. nasa.gov/fap/documents.html, 2006. 

2perry, B., Silva, W., Florance, J., Wieseman, C., Potozky, A., Sanetrik, M., Scott, R., Keller, D., Cole, S., 
and Coulson, D., “Plans and Status of Wind Tunnel Testing Employing an Aeroservoelastic Semispan Model,” 48th 
AIAA/ASME/ASCE/AHS/ASC Structural Dynamics, and Materials Conference, No. AIAA-2007-1770, April 2007. 

^Kopasakis, G., Connolly, J., Paxson, D., and Ma, P., “Volume Dynamics Propulsion System Modeling for Supersonics 
Vehicle Research,” Journal of Turbomaehinery , Vol. 132, No. 4, October 2010, pp. 8. 

‘^Connolly, J., Kopasakis, G., and Lemon, K., “Turbofan Volume Dynamics Model for Investigations of Aero-Propulso- 
Servo-Elastic Effects in a Supersonic Commercial Transport,” 4Sth AIAA/ASME/SAE/ASEE Joint Propulsion Conferenee 
and Exhibit, No. AIAA-2009-4802, 2009. 

^Connolly, J., Kopasakis, G., Paxson, D., Stueber, E., and Woolwine, K., “Nonlinear Dynamic Modeling and Controls 
Development for Supersonic Propulsion System Research,” 4'^th AIAA/ASME/SAE/ASEE Joint Propulsion Conferenee and 
Exhibit, No. AIAA-2011-5653, 2011. 


NAS A/TM— 20 15-21 8479 


15 


^Kopasakis, G., Connolly, J., and Seidel, J., “Propulsion System Dynamic Modeling of the NASA Supersonic Concept 
Vehicle for AeroPropulsoServoElasticity,” 50th AIAA/ASME/SAE/ASEE Propulsion and Energy Eorum, No. AIAA 2014-, 
2014. 

^Chavez, F. and Schmidt, D., “Analytical Aeropropulsive/Aeroelastic Hypersonic - Vehicle Model with Dynamic Analysis,” 
Journal of Guidance, Control, and Dynamics, Vol. 17, No. 6, December 1994, pp. 1308-1319. 

® Nguyen, N., “Integrated Flight Dynamic Modeling of Flexible Aircraft with Inertial Force-Propulsion- Aeroelastic Cou- 
pling,” 44th AIAA Aerospace Sciences Meeting and Exhibit, No. AIAA 2008-194, Reno, Nevada, 7-10 January 2008. 

^Clark, A., Wu, C., Mirmirani, M., and Choi, S., “Development of an Airframe-Propulsion Integrated Generic Hypersonic 
Vehicle Model,” 4^th AIAA Aerospace Sciences Meeting and Exhibit, No. AIAA 2006-218, Reno, Nevada, 9-12 January 2006. 

^^Biedron, R. T., Derlaga, J. M., Gnoffo, P. A., and et. al, “FUN3D Manual 12.4,” Tech. Rep. TM-2014-218179, 2014. 

Anderson, J., “Computational Fluid Dynamics: The Basics with Applications,” McGraw Hill, New York, NY, 1995. 

^^Chang, S.-C., “The Method of Space-Time Conservation Element and Solution Element - A New Approach for Solving 
the Navier-Stokes and Euler Equations,” Tech, rep., NASA TM-111237, 1995. 

^^Connolly, J. and Kopasakis, G., “Loop Shaping Control Design for a Supersonic Propulsion System Model Using QFT 
Specifications and Bounds,” 4^th AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit, No. AIAA-2010-7086, 
2010. 

^^Lytle, J. K., “The Numerical Propulsion System Simulation: A Multidisciplinary Design System for Aerospace Vehicles,” 
Tech, rep., NASA TM 1999-209194, 1999. 

^^Seldner, K., Mihaloew, J., and Blaha, R., “Generalized Simulation Technique for Turbojet Engine System Analysis,” 
Tech, rep., NASA TN-D-6610, 1972. 

Kopasakis, G., “Feedback Control Systems Loop Shaping Design with Practical Considerations,” Tech, rep., NASA TM 
2007-215007, September 2007. 

^^Horowitz, L, “Quantitative Feedback Theory,” IEEE Proceedings , Vol. 129, No. 6, Nov. 1982, pp. 215-226. 


NAS A/TM— 20 15-21 8479 


16 



